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Abstract 

The three invariants of C 1 / 2 are key to expressing this tensor and its inverse as 
a polynomial in C. Simple and symmetric expressions are presented connecting the 
two sets of invariants 12,13} and {11,12,^3} of C and C 1 / 2 , respectively. The 
first result is a bivariate function relating I\,I<i to i\,%2- The functional form of 
i\ is the same as that of 12 when the roles of the C-invariants are reversed. The 
second result expresses the invariants using a single function call. The two sets of 
expressions emphasize symmetries in the relations among these four invariants. 

1 Introduction 

We consider relations among the basic tensors of three dimensional continuum mechanics, 
all defined by the deformation F, 

F = RU = VR, C = F'F, B = FF*. 

U and V are are symmetric and positive definite, and therefore 

U = C 1/2 , V = B 1 / 2 . 

Here we will only consider properties of U and C, but the results apply to V and B. 

Although the square root of a second order positive definite symmetric tensor is unique 
and unambiguous it is not, however, a simple algebraic construct. One way to circumvent 
this problem is to express U as a polynomial in C using the Cayley-Hamilton equation, 

U 3 - ^U 2 + * 2 U - * 3 I = 0. (1) 

Here ii,i2,i% are the invariants of U, 

ii = trU, z 2 = ^(trU) 2 - itrU 2 , z 3 = detU, 

Multiply ([1}) by (U + z'il), and note that the result conta ins te r ms pr oportional to I, U, 
U 2 and U 4 . Replacing the latter two by C and C 2 gives ( Ting! . [l985l ) 



U = (hi 2 - 2a)- 1 (W + (i? - i 2 )C - C 2 ) . (2) 
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Note that — 13 = detail — U) > f Carroll! I2004J ) . The inverse U _1 may be obtained 
by multiplying each side of ([2]) with C" 1 and using the Cayley-Hamilton equation for C 
to eliminate the single remaining C _1 term. The orthogonal rotation tensor follows as 
R = FU _1 , from which o ne can determine kinematic quantities such as the rotation angle 



and the axis of rotation ( Guan-Suol . Il998l ). 

The relation (j2J) for U in terms of C avoids the tensor square root difficulty but 
introduces another: how to express {ii,i2, 23} in terms of C, or more specifically in terms 
of its invariants, 



h = tr C, 



detC. 



While the relation z 3 = is simple, formulas for i\ and %2 are not. But as eq. §Z§ and 
related identities illustrate, the functional relations between the two sets of invariants are 
important for obt aining semi-explicit expressions for stretch and rotation ten sors, and for 
their derivatives ( Hoger and Carlson! 1984a : Steigmann . 120021; ICarroll! 12004 ). 

The first such relations are due to lHoger and Carlson ( 1984bl ) who derived expressions 
for {11,12} by solving a quartic equation. Sawyers! "( 19861 ) subs equently showed that one 
can o btain alternative relations using the standard solutions ( Goddard and Ledniczkyl . 
19971 ) for the cubic equation of the eigenvalue of C. Let Ai,A 2 ,A 3 be the (necessarily 



positive) eigenvalues of U, then the eigenvalues of C are Af, X^, A§, and 



it = Ai + A 2 + A 3 , % 2 ■- 
I\ = A^ + A2 + A3, 1 2 



A1A2 + A2A3 + A3A1, 

■ \\\\ + ^2^3 + ^3^1 > 



v 3' 

h 



A1A2A3, 
= A1A2A3. 



(3a) 
(3b) 



Sawyers' approach is to essentially compute the eigen values of C , take their square roots 
and from these determine the invariants of U by (!3a|) . iJogi ( 2006 ) generalized this scheme 
to tensors of order n > 3. This met hod does not provi d e direc t relati ons between the 
invariants. Although the formulas of iHoger and Carlson! ( 1984b! ) and of I Sawyers! ( 19861 ) 
are explicit, they are not totally satisfactory. In each case the functional forms are com- 
plicated. As we will see, there is no way to avoid this "complexity" since we are dealing 
with roots of cubic and quartic equations. But that is not t he basic issue, rather it is a 
lack of any underlying symmetry or balance in the solutions of lHoger and Carlson! (Il984ri ) 
and of ISawye"rsl ( 19861 ). This makes it difficult to comprehend the formulas, and to place 
them in context. It is all the more unsettling by virtue of the fact that the formulas are 
associated with algebraic systems of deformation tensors, systems that are elegant and 
generally quite transparent. 

The object of this paper is to express {11,12, ^3} in terms of {Ii,l2,Is} in two forms 
that each display the underlying symmetry of the relations. Both forms employ a single 
function, but have slightly different properties. We begin in section [2] with a summary of 
the principal results, followed by a review of the previously known solutions in section [3j 
The new formulas for the invariants of C 1//2 are derived in section HI with some closing 
comments in section [5j 



2 Principal results 

The first result is 
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Theorem 1 The invariants of C 1 ^ 2 



are 



r i/6,/ h h \ . ri/3 w h h \ . r i/2 / i \ 

l l = V/ (7173 .7273)' Z2 = / 3 / (7275' TITS-)' l 3=V, (4) 

where f is a function of two variables, 

f(x, y) = g(x, y) + \/x- g 2 (x, y) + 2/g(x, y), (5a) 

9{x, y) = Q (x + Jx^Zy [(( + Vc^l) l/3 + (C - v/^ 3 !) 173 ] )) ^ (5b) 
27 + 2x 3 - 9xy 

C 2(x 2 - 3y) 3 / 2 ' 1 j 

The function (7 can be expressed in the alternate form 

g(x,y) = )J^(x + 2y / x 2 -3y cos(^ arccosC(x,y))). 

It is clear from Theorem [1] that the following reduced quantities are the important 
variables: 

Ji = TT73' ^ = ^73, (6a) 

^3 ^3 

Ji = 7173, h = 7^, (6b) 

J 3 1 3 

in terms of which the Theorem states 

jl = f(Jl,J2), ji = f(J2,Jl)- 

Alternatively, the sum and difference of reduced invariants may be considered as the 
key parameters, which is evident from: 

Lemma 1 The following relation holds between the invariants of C and C 1 ^ 2 : 

4^=*+* + 2. (7) 
Ji -32 

An immediate consequence of Lemma [U is that we need only determine j\ + j 2 or ji — j 2 
since the other follows directly from ([7j). For instance, we could calculate ji + j 2 = 
/( Ji, J2) + /(^2j -A), but this requires evaluation of / twice, and it does not reveal the 
underlying symmetry of the arguments. The second result is a simpler relation between 
the invariants, one that uses a single call to the function /: 

Theorem 2 The reduced invariants o/C 1//2 and C are connected by 

s J± — J 2 , . 

Jl = 2 + "27TT' (8a) 
32 = 2 " "27TT' (8b) 

s = F(J u J 2 ), (8c) 



where 



by \ tn< , M/3 w 6 + x + y 9 + 5x + 5?/ + x^ 
F(x,y) = (2 + x + y)/ /( (2 + a? + y)2/3 , (g - - + /3 )■ (9) 
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Thus, j\ + j ; 2 = F(Ji,J 2 ) and F is a symmetric function of its arguments, F(x,y) = 
F(y,x). Also, F provides an alternative expression for the function /: 



f(x,y) = 7;F(x,y) + 



2 v 2F(x,y)+A 

This form for / employs the function itself, but evaluated at different arguments. This is 
a property of the nonlinear nature of the function. 



3 The methods of Hoger and Carlson and of Sawyers 

Starting from the identities it may be easily verified that ( Hoger and Carlson! . 1 1984bh 



it-2i 2 = h, i^-2iiz 3 = /2, 4 = h. (10) 

Th e last implies i 3 = iV 2 ■ It r emains to find i\ and i 2 . 

Hoger and Carlsonl ( 1984bl ) eliminated i 2 between eqs. (fTUjb and ffTUl) ?. to obtain a 



quartic equation in ^ which they then solved. The same solution for is obtained more 
directly by starting with the ansatz 

U = \ + P, (ll) 

where A is any one of the triplet {Ai, A 2 , A 3 }. For instance, if A = Ai then p = A 2 + A 3 
and i 2 = Ai(A 2 + A 3 ) + A 2 A 3 is 

* 2 = pA + z 3 /A. (12) 
This holds no matter which value A takes. Substituting from ( fill) into (TlOlh implies 

p 2 = h - X 2 + 2i 3 /\. (13) 

The right member is necessarily positive, and using i 3 = I^ 2 we can therefore express 
p > in terms of I±, J 3 and A. 
In summary, 



ii =A+ yjh -X 2 + 2^/l 3 /X, (14a) 



i 2 = V^a/A + V^iA 2 - A 4 + 2^I 3 X, (14b) 

*3 =V^3, (14C) 

where A is any positive root of the characteristic equation of C, 

A 6 - AA 4 + J 2 A 2 - I 3 = 0. (15) 

For instance, 

a = Q (a + + - (a 2 - 3/ 2 ) 3 ] 1/3 + - v/e - (a 2 - 3/ 2 )3] 1/3 ) ) 1/2 , 

and 

e = i(2/ 1 3 -9/ 1 / 2 + 27/ 3 ). 



Note that we assumed that A in the ansatz (fTTj) and (|T2l) is a root of eq. (fl5j) . but this is 
actually a requirement, as can be seen from eqs. (fTUjb using (TTTj) - (TT3|) . Equations (TTTT) 
and (fT2|) represent a standard method of reducing a quartic t o a cubic equation. 

Equation (I14ap is essentially the same as the first relation^ (jHoger and Carlsonl.ll984bl. 
eq. (5.5 )), although they did not identify the root of the cubic e xplici tly iHoger and Carlson 
f 1984bl ) recommended using ( [101) 9 to obtain i 2 . The relation ( 114bl) is quite different and 
is suggestive of the symmetry underlying the solutions for i\ and i 2 that is evident in 
Theorem [TJ We discuss this further in the next section from a different perspective. 

It is interestin g to comp are this with th e explicit positive solution of ( 1X5]) provided by 
Guan-Suol f 19981 ). based on ISawyersl f 19861 ). Starting with the characteristic equation for 
U, 



A - iiA + i 2 X - h = 0, 



(16) 



combined with eos. (TT0|) 9 a nd f|T0l) q , this becomes a quadratic equation for i 2 . The 
solution is f Guan-Suol Il998l . p. 199) 



\- l {^h + p^h\i + I 2 X>-h). 



This appears to be different than eq. fll4bl) . but they are equivalent when one takes into 
account that A satisfies eo. ([15]) . 

In short. IHoger and Carlson! ( 1984bl ) and lSawyerd ( 19861 ) derived eqs. (I14al) and (114bj) . 
respectively. They did not however, note the symmetry between the formulas, which is 
one of the central themes in this paper: that a single function determines both %\ and i 2 . 
In the next section we complete the proof of Theorem [TJ 



4 An alternative approach 

The three conditions in (flOl) can be combined into a single polynomial identity, 



i 2 z 2 ) 2 + (i lZ - i 3 z 3 ) 2 



hz 4 + I 3 z 6 , 



l + hz 

Using the reduced variables of flU]), this becomes 

(1 - j 2 z 2 ) 2 + (j lZ - z 3 ) 2 = 1 + J lZ 2 + J 2 z 4 + z 6 , \/ze 
Comparing coefficients implies the pair of coupled equations 

jl - 2ji = J2- 

Thus, solutions must be of the form 

jl=f(J\,J2), j2 = f(J2,Jl), 

for some function f(x,y) which satisfies 

f 2 (x,y)-2f(y,x)-x = 0. 



(17a) 
(17b) 



(1* 



1 It should be noted that the second relation in their eq. (5.5) never applies. That is, it can be shown 
that the possible equality cannot occur. 
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This is the fundamental equation for f(x,y). It implies the dual relation 



f 2 (y,x)-2f(x,y)-y = 0. (19) 
Eliminating f(y,x) between eqs. ( 1T81) and ( 1T91) leads to a quartic in / = f(x,y): 

(f 2 - xf - 8f - Ay = 0. (20) 



This is equivalent to the quartic of iHoger and Carlson! ( 1984bl ) but expressed in the 



reduced variables. We have already derived a solution of the quartic in the previous 
section by using the ansatz ( TTTj) based on a root of the cubic equation (fT5l) . Equation 
(I14aj) therefore defines the function /, which can be read off by converting to the reduced 
variables j\,j2, Ji, Ji- It may be easily verified that the function of eq. (jSJ) results. 

But what about the relation (114b|) for z 2 ? It does not seem to convert into the expres- 
sion claimed in Theorem!]], i.e. j 2 = f(J2, Ji)- Rather, using eq. (114bl) and j 2 = /( J2, J\) 
to define / we obtain a different expression for /: 



f(v> x ) = — ( n + V x ~ 9 2 (x,y) + 2/g(x,y)g(x, y). (21) 

9{%i y ) 

This is in fact consistent with the definition of / in Theorem [1] because g{x,y) satisfies 
the normalized version of eq. (1161) . 

g 3 (x, y) - f(x, y)g 2 (x, y) + f(y, x)g{x, y)-l = 0. (22) 

Using this and the expression for f(x,y) in (jSJ), gives eq. ( 1211) . This completes the proof 
of Theorem [TJ 

It is interesting to note from (122]) that l/g(y, x) satisfies the same equations as g(x, y), 

i.e. 

9~ 3 (y, x) - f(x, y)g~ 2 (y, x) + f(y, x)g~ l (y, x) - 1 = 0. 

But this does not mean that g(y,x) equals l/g(x,y), since they can (and do) correspond 
to different roots of the cubic. 

The identity (J7J) of Lemma ([T]) follows from the coupled equations (fTTj) . and the details 
of the proof of Theorem [2] are in the appendix. 



5 Conclusion 

Although the expressions for i\ and 12 involve the roots of the characteristic cubic equation 
of C, it seems that the governing quartic equation ( l20i) is more fundamental. This is the 
equation that defines the functions / and F of Theorems Q] and [2j In fact F is defined by 
/, which is in some ways th e central function involved. It is interesting that the quartic 
equation first considered by IHoger" and Carlson! (jl984bl ) reappears in this manner. 

Which of the expressions for i% and z 2 are actually best in practice? While the expres- 
sions (|SJ) are perhaps the most aesthetically pleasing in form, eq. (JIJ) is probably simpler 
to implement. The final choice is of course left to the reader. 
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Appendix 



A Proof of Theorem [2 



For simplicity of notation, let / and /' denote f(x,y) and f(y,x), respectively. Then the 
coupled equations (fl8"l) and (fT9l) are 

f - V = x, 
f -2f = y. 
Adding and subtracting yields, respectively, 

(/ + /'-i) 2 = i+x +2 /+2/r, 

(/-/ / )(/ + / / + 2)=x-y, 

which in turn imply 

f + f' = F, (A.3a) 

/ " /' = f^f , (A.3b) 

where 

F = 1 + v/l + x + y + 2/f. 

The function F = F(x, y) is clearly a symmetric function of x and i.e., it is unchanged 
if the arguments are switched. Solving eqs. ( 1A.3I) gives 

F x-y_ 

J 2 2F + 4' v ; 

f = — — -^—L. (A.4b) 
17 2 2F + 4 v ; 

Although the formulas (1A.4|) clearly split / into parts that are symmetric and asym- 
metric in the two arguments, they are not explicit since the function F involves the 
product //'. Taking the product of the two expressions in (IA.4I) leads to an equation for 
//'. It is simpler to consider the equation for F, which after some manipulation may be 
reduced to the quartic: 

[F 2 - (6 + x + y)] 2 - 8(2 + x + y)F - 4[(5 + x)(5 + y) - 16] = 0. 

Let 

F = (2 + x + y) 1/3 u, 

then u satisfies 

{u 2 -X) 2 -8m-4F = 0, (A.5) 

where 

6 + x + y 9 + 5x + 5y + xy 

~ (2 + x + y) 2 / 3 ' ~ (2 + x + y) 4 / 3 

Equation (1A. 5[) is the same as the quartic (120!) satisfied by /, but with X and Y instead 
of x and y. Thus, 

u = /(x,y), 

which completes the proof of Theorem [21 
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